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Abstract. We present results of hydrodynamic simulations of star formation triggered by cloud- 
cloud collisions. During the early stages of star formation, low-mass objects form by gravitational 
instabilities in protostellar discs. A number of these low-mass objects are in the sub-stellar 
,-0 ■ mass range, including a few objects of planetary mass. The disc instabilities that lead to the 
^ ' formation of low-mass objects in our simulations are the product of disc-disc interactions and/or 
interactions between the discs and their surrounding gas. 
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^ ' 1. Introduction 

In the classical picture (Shu, Adams & Lizano 1987) stars form from the collapse 
of a single gaseous core in isolation: all stages of the star formation process (i.e. ob- 
i jects of Class 0, I, II, III) as well as planet formation are terminated before the newly- 

formed star/planetary system interacts with its environment (e.g. other stars/planetary 
systems, and/or its surrounding gas). In the more dynamic star formation paradigm of 
gravoturbulent fragmentation (Mac Low & Klessen 2004) interactions during the forma- 
tion process are not only a common phenomenon but also a sufficient condition for the 
formation of binary and higher-order multiple systems, as well as for explaining the shape 
■ of the Initial Mass Function (Jappsen et al. 2005; see also Bate, Bonnell & Bromm 2003, 
, Bate & Bonnell 2005), the evolution of the angular momentum distribution of protostars 
(Jappsen & Klessen 2004), their mass accretion rates (Schmeja & Klessen 2004), the spa- 
tial distribution of young stellar objects (Schmeja & Klessen 2006), etc. The competing 
paradigm of a more isolated star formation process advocated by Krumholz & McKee 
(2005), uses large-scale turbulent support against collapse in a rather isotropic way, giv- 
ing rise to supercritical cores with a variety of masses (sec e.g. Krumholz, McKee & 
Klein 2005 vs. Bonnel & Bate 2006). In contrast, gravoturbulent star formation is based 
on the small-scale density fluctuations produced by supersonic turbulence, with gravity 
then promoting the strongest of these fluctuations, determining the onset of collapse. 

In the context of gravoturbulent fragmentation, Chapman et al. (1992), Turner et al. 
(1995), Whitworth et al. (1995), and Bhattal et al. (1998) have studied star formation 
triggered by cloud-cloud collisions. In essence, the cloud-cloud collision mechanism de- 
scribes the small-scale evolution of the interaction between colliding flows (Vazquez- 
Semadeni et al. 2006, 2007). We note that several authors have studied the interaction 
between large-scale colliding flows in non-self-gravitating media, discussing the impor- 
tance of such interactions for the formation of molecular-cloud structure and turbulence 
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(e.g. Hunter et al. 1986, Walder & Folini 1998, Heitsch et al. 2005, Hennebelle & Audit 
2007). These authors discuss the formation at best of low-mass structures such as dense 
cloud cores, whereas in our self-gravitating simulations we use such dense cores as the 
initial conditions in order to follow star formation. 

In this paper we present the results of recent high-resolution hydrodynamic simulations 
of collisions between low-mass clumps (Kitsionas & Whitworth 2007). Due to the high 
numerical resolution of our simulations, we can resolve, for the first time in simulations 
of cloud-cloud collisions, the formation of sub-stellar objects, some of which have a mass 
that lies close to the boundary between brown dwarves and planetary-mass objects, i.e. 
~ 15Mj. These planetary- mass objects form by gravitational instabilities in the disc of 
their parent protostar, i.e. by the same mechanism that is responsible for the formation 
of low-mass stellar and/or brown dwarf companions to the central protostar of such 
discs. Moreover, the disc instabilities that lead to the formation of low-mass objects in 
our simulations are the product of disc-disc interactions and/or interactions between 
the discs and their surrounding gas. Furthermore, we deal here with disc instabilities of 
young, massive discs still in the process of formation, i.e. at the early stages of the disc 
life. We therefore caution the reader not to confuse the instabilities reported here with 
the extensive work on gravitational instabilities and related phenomena in thin, or in 
general evolved discs that have been recently reviewed by Durisen et al. (2007). As for 
the core accretion model of planet formation, we refer the reader to a number of papers 
in this volume dealing with the different stages of solid core formation, gas accretion and 
migration in protoplanetary discs. 

In section 2, we give a brief description of the model and the numerical method we 
use. In section 3, we present our results. In section 4, we discuss our findings with respect 
to models of gas giant planet formation due to gravitational instabilities in protostcllar 
discs. 

2. Our model and numerical method 

In Kitsionas & Whitworth (2007), we investigate, by means of numerical simulations, 
the phenomenology of star formation triggered by low-velocity collisions between low- 
mass molecular clumps. The simulations are performed using a smoothed particle hydro- 
dynamics (SPH) code which satisfies the Jeans condition by invoking on-the-fly Particle 
Splitting (Kitsionas & Whitworth 2002). 

Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebcrt, 
spheres. Collisions are characterised by Mq (clump mass), b (offset parameter, i.e. ratio of 
impact parameter to clump radius), and M. (Mach Number, i.e. ratio of collision velocity 
to effective post-shock sound speed). The gas follows a barotropic equation of state, 
which is intended to capture (i) the scaling of prc-collision internal velocity dispersion 
with clump mass (Larson 1981), (ii) post-shock radiative cooling, and (iii) adiabatic 
heating in optically thick protostellar fragments (with exponent 7 ~ 5/3). The equation 
of state we use is given by Eq. 1 of Kitsionas & Whitworth (2007) and is graphically 
illustrated in the temperature-density plane by their Fig. 1. In short, the initial effective 
isothermal sound speed of the gas, i.e. when non-thermal pressure due to turbulence is 
included, follows the Larson (1981) relations in order to model the observed variation 
of the internal velocity dispersion with clump mass. As long as the collisions begin, the 
effective sound speed is reduced according to a P oc p 1 / 3 relation until it reaches a value 
equivalent to the temperature of 10 K. Then the gas remains isothermal with increasing 
density up to the point where it becomes opaque to its own cooling radiation and heats 
up as P cx p 5 / 3 . 
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Figure 1. Column density plots for the collision with Mo = WMq, b = 0.2, M — 10. Left Panel. 
Initial conditions viewed along the z-axis; Ax = 0.92 pc, Ay = 0.56 pc; the grey-scale is logarith- 
mic, in units of g cm -2 , with sixteen equal intervals from l.OOx 10~ 3 gem -2 to 2.69 x 10~ 2 g cm -2 . 
Right Panel, t — 0.396 Myr; viewed along the z-axis; Ax = 0.016 pc, Ay — 0.014 pc; sixteen-in- 
terval logarithmic grey-scale, in units of g cm -2 , from 2.14 x 10" 1 gem -2 to 3.55 x 10 3 gem -2 . 



From the variety of models we have investigated, we present here the collision between 
two Mo = 10M© clumps, each moving with velocity (along the x-axis) w c i ump = 1 km 
s -1 (corresponding to a Mach number M. = 10, i.e. assuming a post-shock sound speed 
of 0.2 km s _1 ), and with an offset parameter b = 0.2 taken along the y-axis. The initial 
conditions are shown in the left panel of Fig. [T] 

We evolve the simulations using an SPH code. This is a Lagrangian method for 3D 
hydrodynamics. In SPH the sampling points (particles) are able to move with the fluid 
and their properties are distributed in space through a smoothing function that allows 
all hydrodynamic quantities to be continuous (Monaghan 1992). Our self-gravitating 
SPH code (Turner et al. 1995) uses tree-code gravity (Barnes & Hut 1986) to reduce the 
computational cost for the calculation of gravitational forces, and a second order Range- 
Kutta time integrator with multiple particle timesteps. It also uses Particle Splitting 
(Kitsionas & Whitworth 2002), which allows the on-the-fly increase of numerical resolu- 
tion (in terms of splitting the SPH particles) when and where this becomes necessary, i.e. 
every time that violation of the Jeans condition (Bate & Burkert 1997, Truelove et al. 
1997) becomes imminent. The benefit of the use of Particle Splitting is twofold. Firstly, 
gas fragmentation remains resolved at all times. Secondly, the simulations evolve with 
minimum computational cost, as high numerical resolution is employed only when re- 
quired. Our Particle Splitting method assigns velocities to new particles by interpolating 
over the velocity field of their parent particle population, thus achieving conservation of 
energy and momenta with very high accuracy. 

Our code includes sink particles (Bate, Bonncll & Price 1995). These are collisionless 
(star) particles that replace dense, collapsed gaseous regions in order to allow the sim- 
ulations to evolve further in time. With sink particles, the simulations avoid using tiny 
timesteps to follow the evolution internal to the sink but at the same time all information 
on scales smaller than the size of the sink particle becomes invisible to us. In this work 
we use sink particles with radii ~ 20 AU. A summary of all features of our code can be 
found in Kitsionas & Whitworth (2007). 
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Figure 2. Column density plots for the collision with Mo = WMq, b — 0.2, M = 10, with 
sink particles overlayed as solid circles (note that the symbol used for the sinks is larger than 
the actual sink radius at the scale of these plots). View along the y-axis; Ax = Az — 0.02 pc 
in both panels. Left Panel, t ~ 0.393 Myr, i.e. very close to the time of the right panel of 
Fig- HI sixteen-interval logarithmic grey-scale, in units of g cm" 2 , from 3.31 x 10" 1 gcm -2 to 
5.25 x 10 2 gcm~ 2 . Right Panel, t ~ 0.402 Myr; sixteen-interval logarithmic grey-scale, in units 
of g cm" 2 , from 2.88 x lO^gcm" 2 to 2.04 x 10 2 gcm~ 2 . 



3. Results 

A shock-compressed layer forms at the collision interface as the collision proceeds. 
The right panel of Fig. Q] offers a zoom-in view of this layer at t ~ 0.396 Myr. The layer 
extends along the z-axis. The x-z projection of the shock-compressed layer at this time is 
shown in the left panel of Fig. [2j The layer has broken up into filaments and the filaments 
have fragmented into protostars surrounded by circumstcllar discs. Sink particles have 
replaced the central protostars of these discs. A second protostar forms on the upper part 
of the filament that extends almost vertically in this panel. At t ~ 0.402 Myr (right panel 
of Fig. [J) the two protostars in this filament (represented by sink particles surrounded 
by circumstellar discs) have moved along the filament towards each other. Due to the 
tumbling nature of the filament, they have a close encounter in the course of which 
they capture each other into a binary. During this interaction, their discs merge to form 
a circumbinary disc (left panel of Fig. [3]). This disc-disc interaction leads also to the 
fragmentation of the circumbinary disc (see the additional sink particles in this panel). 
At later times, the circumbinary disc fragments only when lumps of material from the 
filament fall on to the disc (e.g. see the lumps on the right panel of Fig. [3]). 

Each of the low-mass fragments, which form during the disc-disc interaction and/or 
during the episodes of lumpy accretion from the filament, gets ejected from the system 
after a close encounter with the binary. The final state of the system (Fig. [4]) includes 
a large circumbinary disc hosting a binary, with components having mass ~ l.5M@ 
(each) and separation < 40 AU, as well as a number of low-mass objects that have been 
ejected from the disc with moderate velocities of order a few km s _1 . Some of these 
low-mass objects are sub-stellar or even of planetary mass. After their ejection from the 
circumbinary disc they have stopped accreting more mass. 
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Figure 3. Column density plots for the collision with Mo = IOMq, b = 0.2, M = 10, with sink 
particles overlayed as solid circles (note that the symbol used for the sinks is larger than the 
actual sink radius at the scale of these plots). View along the y-axis; Aa; = Az = 0.02 pc in both 
panels (note that these panels have been shifted by 0.008 pc along the z-axis with respect to the 
panels of Fig. [2}. Left Panel, t ~ 0.417 Myr; sixteen-interval logarithmic grey-scale, in units of g 
cm" 2 , from 2.45 x 10" 1 gem -2 to 9.12 x 10 1 gem" 2 . Right Panel, t ~ 0.435 Myr; sixteen-interval 
logarithmic grey-scale, in units of g cm" 2 , from 3.09 x 10" 1 gem" 2 to 9.55 x 10 2 gem" 2 . 



4. Discussion 

Based on the results presented here, we can claim with confidence that dynamical 
interactions during the early stages of star formation provide a mechanism for the for- 
mation of low-mass objects. Disc-disc interactions and/or interactions between the discs 
and their gaseous environment (e.g. episodes of lumpy accretion from a filament) can lead 
to gravitational instabilities in protostellar discs and the formation of low-mass objects. 
Some of these objects may subsequently be ejected from the disc they formed in. As a 
result of this ejection, their accretion rates will drop significantly and they will end-up 
in the sub-stellar or even planetary-mass regime. 

It is important to note that the gas in our discs is rather warm, as no cooling mechanism 
for discs is included in our code and adiabatic gas heating has already switched on at the 
gas densities of our discs (see section 2 for a brief description of our equation of state) . 
Therefore, the gravitational instabilities we report here are due to an increase in the disc 
surface density (caused by a dynamical interaction) and not because of the decrease of 
the gas sound speed in the disc. 

We note that the discs in our simulations contain only a few thousand SPH parti- 
cles and evolve for a small number of orbital periods. Thus, due to excess numerical 
dissipation and limited orbital evolution, we are not yet in a position to investigate in 
detail the angular momentum evolution of gas and fragments in our discs. Nevertheless, 
fragmentation here is mostly externally triggered and therefore we are confident that 
it is not due to spurious angular momentum transport. Moreover, since the discs form 
self-consistently in a cluster environment, it is inevitable that during the initial stages of 
their formation they will contain a small but increasing number of particles. This makes 
impossible to comment on the effect of numerical dissipation on early disc evolution. 
Convergence studies with respect to the disc numerical dissipation in our simulations are 
currently underway. 

Gas cooling is an important ingredient in simulations of isolated protostellar discs 
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Figure 4. Column density plots for the collision with Mo = IOMq, b = 0.2, M. = 10, with sink 
particles overlayed as solid circles (note that the symbol used for the sinks is larger than the 
actual sink radius at the scale of these plots). The final state of the simulation (t ~ 0.440 Myr); 
view along the y-axis; Aa; = Az = 0.02 pc; sixteen-interval logarithmic grey-scale, in units of g 
cm~ 2 , from 3.39 x lO^gcm -2 to 6.61 x 10 2 gcm~ 2 . 



(e.g. see Rice et al. 2003) even in the case that this cooling is not sufficient to produce 
gravitational instabilities and fragmentation of the disc (see the discussion in Durisen 
et al. 2007). 

Stamatellos, Hubber & Whitworth (2007) have recently reported gravitational instabil- 
ities that lead to the formation of low-mass objects in an isolated disc using an SPH code 
with on-line (approximate) radiation transport. Their radiation hydrodynamics method 
is effective in calculating disc cooling close to the disc midplane but not further out in 
the disc "atmosphere" . We are currently developing an improved algorithm that will be 
able to identify the scale-height at which each gas particle lies and thus assign the cooling 
efficiency of the particle accordingly. 

Generalising our finding on the importance of dynamical interactions during the early 
stages of star formation, we assert that such interactions may be important for the 
formation of low-mass objects also at later stages of protostellar evolution, as long as a 
gas disc is still present. The outcome of such interactions need not be free-floating objects 
in all cases, as in less "chaotic" situations (e.g. in circumstellar discs around single stars) 
the low-mass objects formed may remain bound to the system they formed in. Moreover, 
interactions of already formed planetary systems with incoming stars or lower mass free- 
floating objects may disturb the planetary systems, but this is beyond the scope of this 
paper (see several papers in this volume on planetary system dynamics) . 

Taking our results at face value, one can claim that the majority of the low-mass objects 
forming in such simulations will end up unbound from the system of their birth, i.e. a 
number of free-floating objects are expected to form as a by-product of star formation. 

The formation of free-floating planetary-mass objects is of particular interest. Such 
objects can be observed mainly through their gravitational signature. In the field, this 
signature could be detected only through microlensing observations. Zinnecker (2001) 
has calculated that a lensing planet can produce a microlensing event with duration of a 
few days. He has also estimated that one such event is expected out of 3 x 10 s observed 
lightcurves. He has further advocated that campaigns for the observation of free-floating 
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planetary-mass objects through microlensing can become possible when the VST and 
VISTA facilities start their operation as well as with the JWST. 

The event probability and duration that Zinnecker (2001) has estimated are based on 
simple calculations involving the known frequency of binary stars with members hosting 
planetary systems and the probability for dynamical interactions between the members 
of such complicated systems. It would be interesting to calculate the corresponding num- 
bers based on our models. However, we need conduct larger simulations of gravoturbulent 
star formation (see e.g. Klessen & Burkert 2000, 2001, Klessen 2001a, b, Jappsen et al. 
2005) in order to be able to derive more reliable statistics on the number of free-floating 
planetary-mass objects as well as the overall mass distribution of objects forming in the 
calculations. Such simulations are currently underway. Moreover, by using multiple levels 
of Particle Splitting we are now able to employ even higher numerical resolution for our 
discs as well as use smaller sink particles. We are thus in a position to investigate the 
evolution of individual discs within a realistic cluster environment at a resolution compa- 
rable to that of simulations of isolated discs. This way we can now address issues such as 
the amount of numerical dissipation and its effect on the angular momentum evolution 
of our discs. 
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